Unusually strong coherent response from grain-boundary Josephson network in 

polycrystalline Pr x Yi_ x Ba 2 Cu 3 0-7^s 
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By applying a highly sensitive homemade AC susceptibility technique to Pr x Yi- x Ba2Cus07-s 
polycrystals (with x = 0.0, 0.1 and 0.3), we observed very sharp Fraunhofer type low-field periodic 
oscillations of the real part of the AC susceptibility which were attributed to Josephson vortices 
penetrating intergranular regions of grain-boundary Josephson network in our samples. Assuming 
the Lorentz type distribution of single-junction contact areas, we were able to successfully fit the 
experimental data. 

PACS numbers: 74.25.Ha; 74.50.+r; 74.62.Dh 
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INTRODUCTION 

It is well-known that important for large-scale ap- 
plications properties of any realistic device based on 
Josephson effects require a very coherent response from 
many Josephson contacts comprising such a device (see, 
e -g-> 0) 0j S IH IS H, 0i H] an d further references therein). 
Usually, due to inevitable distribution of critical currents 
and sizes of the individual junctions, a grain-boundary 
induced Josephson network in polycrystalline materials 
manifests itself in a rather incoherent way, making it vir- 
tually impossible for applications. That is why, ordered 
artificially prepared (hence more costly) Josephson junc- 
tion arrays (JJAs) are used instead to achieve the ex- 
pected performance 0, [l(J EH, 0, H Q ■ 

In this Letter we report on unusually strong coher- 
ent response of grain-boundary Josephson network in our 
polycrystalline Pr x Y\^ x Ba2Cu^Oi (PYBCO) samples 
which manifest itself through a clear Fraunhofer type 
magnetic field dependence of the measured AC suscepti- 
bility (more typical for ordered JJAs). 



RESULTS 

High quality PYBCO bulk polycrystalline samples 
have been prepared by following a chemical route based 
on the sol-gel method [l5| . The phase purity and the 
structural characteristics of our samples were confirmed 
by both scanning electron microscopy (SEM) and x-ray 
diffraction (XRD) along with the standard Rietveld anal- 
ysis. The analysis of the XRD data (Fig.l) reveals that 
no secondary phases are present in our samples and that 
the peaks correspond to the orthorhombic structure with 
Y Ba2Cuz07-s (YBCO) stoichiometric phase. The on- 
set temperatures Tc (shown in Fig. 3 for all studied sam- 
ples) were independently confirmed via the resistivity, 
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FIG. 1: XRD patterns of Pr x Yi— x Ba2Cu30r-s samples. 



magnetization and AC susceptibility data, and well cor- 
relate with the values reported in the literature for poly- 
crystalline samples with a similar composition 

EE El- 
Fig. 2 shows the SEM scan of grain-boundary morphology 
in the undoped YBCO sample (with grains of different 
shape and average size of the order of lpm). 

AC measurements were made by using a high- 
sensitivity homemade susceptometer based on the screen- 




FIG. 2: SEM scan photography of YBa 2 Cu30 7 - 5 . 



ing method and operating in the reflection configura- 
tion [^, [3] ■ The complex response Xac — X + *X was 
measured as a function of the AC field h ac (t) = hcos(u>t) 
(applied normally to the sample's surface with the am- 
plitude O.OlOe < h < 50Oe and frequency 1kHz < u> < 
30kHz) taken at fixed temperature. The field depen- 
dence of the normalized real part of AC susceptibility 
A\ (/i) = x {h) ~ X (0) f° r different temperatures and 
Pr content is shown in Fig. 3. As it is evident from 
this picture, there are two distinctive regions. Namely, 
above h = 150e the curves exhibit almost linear depen- 
dence which can be attributed to establishment of the 
well-documented [la, |2CJ, l2l| Bean type critical state 
regime with l+47r% B (/i) = 2h/JcD where Jc is the field- 
independent critical current density and D the sample's 
thickness. On the other hand, below h = 150e, prac- 
tically temperature-independent periodic oscillations are 
clearly seen. In what follows, we shall focus on explana- 
tion of this interesting phenomenon. 



DISCUSSION 

To understand the observed behavior of A% (h), it 
is quite reasonable to assume that the low-field AC re- 
sponse in our samples is related to a Josephson network 
mediated intergranular contribution [i^, HH, Xj(h). 
Notice, first of all, that the field oscillations manifest 
themselves in a rather narrow region between two char- 
acteristic Josephson fields: the lower critical field h J cl = 
$o/27rAj and the upper critical field h J c2 = ^q/2-kXjCL. 
Here, Xj = \J ^o/2wfioJcd is the Josephson penetration 
depth (which also defines the size of a Josephson vortex) , 
and d = 2\l + I ~ 2Al is the width of a single contact 
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FIG. 3: The magnetic field dependence of the normalized real 
part of AC susceptibility at different temperatures for samples 
with different Pr content x. Inset: the best fits (solid line) of 
the low-field region according to Eqs.(2)-(4). 
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with Xl being the London penetration depth and I the 
thickness of the insulating layer. Based on the above crit- 
ical fields, the flux penetration scenario can be described 
as follows. For h < h^i, we have analog of the Meiss- 
ner state for intergranular region (total screening of the 
applied field). When h > h%±, field starts to penetrate 
into the intergranular region in the form of Josephson 
vortices (hence, h J cl is analog of the Abrikosov lower crit- 
ical field h£y = §q/2-k\ 2 l ). Since (unlike Abrikosov vor- 
tices) Josephson vortices are coreless, they will nucleate 
in the contact area 5 = Ld (where L is the length of the 
contact) until their distribution becomes homogeneous. 
This process takes place for field region h J c2 > h > h J cl 
and manifests itself in the Fraunhofer type dependence 
of the Josephson current I.j(h) = 7/(0) (sin///) with 
/ = h/hj. The structure of the pattern is governed by 
the characteristic Josephson field hj — 2irS/$o which 
indicates how many fluxons penetrate the contact area 
5. At h ~ h'l 2 , the Fraunhofer pattern practically dis- 
appears. More precisely, Ij(h&) <C lj(0) (hence, h J c2 
indeed plays a role of the upper critical field for Joseph- 
son vortices). For h£y > h > fyf 2 , we have a conventional 
Meissner state for the intragranular region (applied field 
still can not penetrate inside grains). The formation of 
the intragranular Abrikosov mixed state starts only for 
h > hf v 

In order to clarify the origin of the discussed here ef- 
fects, let us estimate the values of the above critical fields 
(and the corresponding depths). By relating the lower 
Josephson field h J cl = $o/27rA 2 with the beginning of 
the observed oscillations (which is h J cl ~ 3.50e for pure 
YBCO sample), we obtain A,/ ~ lfim for the size of the 
Josephson vortex. On the other hand, by relating the 
upper Josephson field h J c2 = <& /27rAjd with the end of 
oscillating behavior (which is h J c2 ~ 150e for the same 
sample), we get Xl — lOOnm as a reasonable estimate 
of the London penetration depth in this material [la |. 
To account for the observed evolution of the Josephson 
fields with Pr concentration x, we recall [l7| that in ad- 
dition to the critical temperatures Tc{x) c± Tc(0)(l — x), 
both the London depth Xl(x) and the critical current 
density Jc(x) decrease upon doping. Namely, assuming 
that X L (x) ~ A L (0)(l-x) and J c {x) ~ J C (0)(1 - a;), we 
obtain Xj(x) ~ Aj(0)/(1 - x), h J cl {x) = <& /2irX 2 j(x) ~ 
h J cl (Q){l - x) 2 and h J c2 {x) = $„/ 'AitXj{x)X l {x) ~ h J c2 (0) 
for doping induced variation of the Josephson depth, 
lower and upper Josephson fields, respectively. Notice 
first of all that (in accord with our observations) the 
upper field does not change with Pr and has the same 
value (h J c2 ~ 150e) for all three samples. As for the 
lower field, according to the above simplified expres- 
sions, we have h J cl {x = 0) ~ 3.50e, h J cl (x = 0.1) ~ 
0.8/1^(2 = 0) ~ 2.80e, and h J cl {x = 0.3) ~ Q.4Qh J cl {x = 
0) ~ 1.750e for x — 0, x = 0.1, and x = 0.3, respec- 
tively. All these estimates are in good agreement with 
the observed onset of field oscillations (see Fig. 3). Be- 



sides, the above cc-dependence of the London penetration 
depth results in the following evolution of the contact 
area S(x) ~ 5(0)(1 - x) where 5(0) ~ 2LX L (0). Now 
we can also estimate the value of this area for each of 
three samples by relating the number of trapped fluxons 
n(x) = 2irhS(x) /$o with the number of the observed os- 
cillation minima (seen in Fig. 3) for the applied field span 
(lying between the lower and upper critical fields) h ~ 
h^ 2 (x) — hci{x). Namely, using n(0) = 4, n(0.1) = 3, and 
n(0.3) = 3, we obtain 5(0) ~ l^m 2 , 5(0.1) ~ 0.9fim 2 , 
and 5(0.3) ~ 0.7 fim 2 for the estimates of contact areas in 
our three samples, which remarkably correlate with the 
above assumed doping dependence of S(x). Moreover, 
the doping dependence of the London penetration depth 
Xl{x) also controls the evolution of the lower Abrikosov 
field, h£_(x) = fc£(0)/(l - a;) 2 leading to the follow- 
ing estimates: h^(0) ~ 0.03T, ^(0.1) ~ 0.04T, and 
/i^\(0.3) ~ 0.06T. Therefore, given a markedly differ- 
ent values for Josephson and Abrikosov critical fields, we 
can safely assume that the discussed here phenomenon is 
strictly related to the Josephson physics. A pronounced 
Fraunhofer type form of the observed curves suggests a 
rather strong coherent response from many Josephson 
junctions comprising the grain-boundary network (de- 
spite some distribution in sizes of the individual junctions 
seen in Fig. 2). To describe the observed phenomenon, we 
assume that intergranular contribution Xj{h) is related 
to AC field h ac (t) = hcos(uit) induced modulation of the 
Josephson current hj{t) = (2n /<&o)Jij sin0y(t) (where 
Jij is the Josephson energy) circulating in a closed pla- 
quette (cluster) with a random distribution over contact 
areas Sij. Each such cluster involves adjacent supercon- 



ducting grains i = 1,2,. 
tive phase difference 0u 
tergranular barriers [^, [23J, |24j, [25 



,N and j = 
2t 



■ i + 1 with an effec- 
h ac {t)/$o across in- 
J . In turn, due to the 
Ampere's law, this circulating current hjit) produces a 
net magnetic moment [H, El] fi(t) = Iij(t)Sij, lead- 
ing us to TL{t) — f° r the total Hamiltonian 
describing the flux dynamics of a single plaquette with 

2n 

Hij(t) = Jij[l-cosdij(t)]- — JijsinOij(t)hac(t)Sij (1) 
where the second term is a Zeeman contribution 

fl(t)hac(t). 

To obtain the experimentally observed intergranular 
contribution to the AC response, we assume (for simplic- 
ity) a Lorentz type distribution of the contact areas Sij 
(around their mean values Sq with the width A) of the 
form: 



F(Sij) 



1 



7T (Sij - S ) 2 + A 2 



(2) 



keeping in mind that the Josephson energy Jij also 
depends on the contact area Sij via a distance be- 
tween grains r^. Namely, according to the conven- 
tional description of granular superconductors, Jij = 



4 



J(0) exp(— Tij/2ro) where ro is of the order of an av- 
erage grain size (radius). By using some geometri- 
cal arguments, it can be demonstrated that SV, ~ 
i r ij/ r o] ) So which results in the following explicit depen- 
dence of the Josephson energy on contact area, = 
J(0) exp(— ^/Sij/ASo) . Notice that this way we do not in- 
troduce any new fitting parameters (apart from the above 
mentioned So and A). 

Thus, the expression for the observed intergranular 
contribution to the AC susceptibility finally reads 

X (h) =< x'j(h) > s = £ I' 3 ™ dS t] F{S l3 ) Xt] {h) (3) 
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where 



2nVj 



d{wt) 



dhl c (t) 



(4) 



= Xij(o)[Mfij) - Ji(fij)h 



Here, x y (0) = (2-rrSij /$ ) 2 Jij /V (V is the properly de- 
fined volume), J n are the Bessel functions, and fij = 
h/h^j, with hfj = $ /27r5ij being a characteristic 
Josephson field (which is eventually responsible for the 
structure of the Fraunhofer pattern). 

The best fits of the experimental data for the low- 
field region based on Eqs.(2)-(4) along with the values 
of two fitting parameters, So and A, are shown as in- 
sets in Fig. 3 (we assume S m — 2So). In particular, 
for undopcd YBCO sample, we found So(0) = 1.12/im 2 
for the mean value of the contact area (which corre- 
sponds to the characteristic Josephson field hj(0) = 
$o/27rSo(0) ^ 2.80e) and A(0) = 0.13^to 2 for the con- 
tact area width distribution. The contact areas in the 
doped samples are found to be best fitted by the fol- 
lowing set of parameters: 5o(0.1) = 0.95^to 2 (equiv- 
alent to MOT) = $ /27t5o(0.1) ~ 3.3<3e), A(0.1) = 
0.16/wn 2 , So(0.3) = 0.81^to 2 (equivalent to hj(0.3) = 
$o/2ttSo(0.3) ~ 3.9<3e), and A(0.3) = 0.22 M m 2 . Notice 
that decreasing of So(x) with x closely follows the earlier 
suggested doping dependence Sq{x) ~ So(0)(l — x), while 
the opposite behavior of the widths A (a;) (increasing in 
doped samples) most likely reflects a random accumula- 
tion of Pr on grain boundaries, leading to a more broad 
distribution of the contact areas. 

In summary, by using a highly sensitive homemade 
AC magnetic susceptibility technique, the magnetic 
flux penetration has been measured in high-quality 
Pr x Y 1 ^ x Ba2Cu2,0'j-s polycrystals as a function of AC 
magnetic field for different temperatures and Pr doping. 
In addition to the conventional critical state behavior at 
higher fields, a clear manifestation of coherent intergran- 
ular response from Josephson vortices seen as a periodic 
Fraunhofer type dependence of the real part of AC sus- 
ceptibility was observed at low magnetic fields. 
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